Integrable versus Non-Integrable Spin Chain Impurity Models 



Erik S. S0rensen( a ), Sebastian Eggert^, and Ian Affleck^' 6 ) 
^Department of Physics and ^ Canadian Institute for Advanced Research 
University of British Columbia, Vancouver, BC, V6T 1Z1, Canada 
(February 1, 2008 UBCTP-93-007) 

Abstract 

Recent renormalization group studies of impurities in spin-1/2 chains appear 
to be inconsistent with Bethe ansatz results for a special integrable model. 
We study this system in more detail around the integrable point in param- 
eter space and argue that this integrable impurity model corresponds to a 
non- generic multi-critical point. Using previous results on impurities in half- 
integer spin chains, a consistent renormalization group flow and phase diagram 
is proposed. 
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I. INTRODUCTION 



Recently there has been considerable interest in various quantum impurity problems.0~i 
These can generally be formulated as one-dimensional Luttinger liquids interacting with 
local defects of various kinds. In general it is expected that such quantum impurity models 
renormalize to critical points which correspond to conformally invariant boundary condi- 
tions. The quantum impurity is screened and/or decouples; it does not appear in the fixed 
point Hamiltonian although a remnant effective impurity, decoupled from the continuum 
degrees of freedom, may be left behind. 

A particularly simple example is a single impurity in a spin S = 1/2 Heisenberg antifer- 
romagnetic chain. Two of the present authors analyzed a large class of models of this type 
using analytic renormalization group (RG) arguments and numerical finite-size analysisi. 
We concluded that the only stable critical points correspond to a completely unperturbed 
chain or else a chain with a break at the impurity location. Taking the initial boundary 
conditions to be periodic, we refer to these two critical points as the periodic and open chain, 
respectively. In the open case, but not the periodic, a remnant impurity spin may also be 
present, as we shall review in section |T[ It was recently drawn to our attention that a Bethe 
ansatz integrable impurity model of this type was solved several years ago. Its low-energy 
behaviour corresponds to a conformally invariant boundary condition, but not to one of the 
stable critical points mentioned above. The purpose of the present article is to resolve this 
apparent contradiction. 

The integrable impurity model involves a single spin-S* impurity which is coupled sym- 
metrically to two neighbouring sites on the chain. The Hamiltonian, found by Andrei and 
Johannessonl is 

i=l 

+ ^(^Tl) 2 ^ -S + a L -S+ 1{B X -S,B L -S}- S(S + l)B x ■ B L \. (1.1) 

where the aVs are Pauli matrices and {, } denotes the anti commutator. In the following 
we shall refer to Eq. ( |1 . 1|) as the integrable impurity model. The equivalent problem of a 
spin-1 chain coupled to a spin-S" impurity was solved by Lee and Schlottmanni, and later 
generalized to a spin-S" chain coupled to a spin-S impurityi. For S > S' it was found 
that in the thermodynamic limit the system behaves like a spin-S" chain with one extra site 
and a decoupled spin S e s of size S — S'. The S' > 1 models already exhibit non-generic 
behaviour before adding the impurity. In particular these integrable periodic models00, 
without impurities, do not exhibit the Haldane gap for integer S'. For a discussion see 
Ref. We shall not consider them further. 

A peculiar feature of the integrable impurity model with S' = 1/2 and S > 1 is that 
the effective, partially screened impurity at the critical point has spin S c s = S — lfM, 
despite the fact that the impurity couples with equal strength to two spin-l/2's. This 
seems contradictory since, if we assume that the critical point corresponds to an infinite 
antiferromagnetic coupling then we would obtain S c g — S — 1. Furthermore, our RG analysis 
indicates that in any event if S e R ^ a stable critical point must correspond to the open 
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chain. Our further analysis of the integrable impurity model with 5 = 1, discussed in Section 
III, indicates that the critical point corresponds to the periodic chain with S c q = 5 — 1/2 = 
1/2. It is as if the impurity spin "splits in half", donating an extra 5 = 1/2 spin to the 
periodic chain and leaving behind a decoupled s = 1/2 impurity. [See Fig. [[[] 

We argue below that this corresponds to an unstable critical point which is peculiar to this 
Hamiltonian. Generic Hamiltonians renormalize to the stable fixed points mentioned above. 
Thus it appears that the conditions for integrability somehow "fine-tune" the impurity-model 
Hamiltonian so that it corresponds to an unstable fixed point. The same phenomenon was 
found earlier for integrable periodic chains of spin 5" > l£3. 

In the next section we briefly review our RG analysis which shows that the integrable 
impurity model cannot correspond to any of the known stable fixed points. We then conjec- 
ture that it corresponds to the particular unstable fixed point mentioned above. In Section 
III we analyze this unstable fixed point. In particular, we find that the RG flow to this fixed 
point is governed by two marginally irrelevant operators which lead to finite-size corrections 
which only go away with the inverse logarithm of the chain length. Fortunately, we are able 
to calculate energy eigenvalues for chain lengths up to 5,000 using the Bethe ansatz. This 
enables us to analyze in detail the logarithmic behaviour and show convincingly that our 
conjecture is correct. In Section IV we analyze the effect of perturbing the couplings to 
the impurity. We conjecture a general RG flow and phase diagram and attempt to test it 
numerically. The non-integrability limits the maximum chain length to about 20. Because 
of the presence of two marginal operators it is difficult to draw definitive conclusions from 
these calculations but they seem to be consistent with our conjecture that the integrable 
impurity model corresponds to an unstable fixed point. 



II. REVIEW 

The continuum limit of the 5=1/2 Heisenberg antiferromagnet, Si = (l/2)<Tj, can be 
written in terms of a free boson, with a particular value of the "compactification radius" (or 
Lagrangian normalization) which enforces the SU(2) symmetry. (For a review see Ref. ||) 
This model is equivalent to the SU(2) Wess-Zumino-Witten model with Kac-Moody cen- 
tral charge, k — 1. The uniform and staggered magnetization correspond to two different 
operators. Thus the spin operators, Si = (l/2)a», becomes: 

Si « (J L + J R ) + (-Inconstant • tvha. (2.1) 

Here Jl and Jr are the left and right-moving spin densities or currents and h is an SU(2) 
matrix field. The current operators have scaling dimension x = 1 while h has x = 1/2. 
Using the operator product expansion, it can be shown that, 

5 4 • S l+1 « constant (-l) J tr/i. (2.2) 

We now review the effect of local perturbations upon the open and periodic chain fixed 
points. 

The various types of local perturbations of the periodic chain corresponding to quantum 
impurities can all be expressed in terms of J = Jl + Jr and h, the former operator being 
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marginal and the latter relevant. In fact, most perturbations generate relevant operators, 
the only exception being perturbations symmetric under site-parity Vs, ie. reflection about 
a site, which do not involve an external spin. In this case the parity symmetry ensures that 
all terms involving trh cancel. A perturbation which is not invariant under Vs will generate 
trh under a renormalization group transformation. Let us now consider perturbations that 
do involve an external spin. A Vs invariant coupling, i.e. ff{ ■ S, to an external spin S 
generates the relevant operator trha ■ S, as can be seen from Eq. Q2.1| ). (Note that tvha ■ S 
is a relevant operator since a decoupled impurity has zero scaling dimension.) For equal 
Heisenberg coupling of two neighbouring chain-spins to the impurity, symmetric under link- 
parity, i.e. d\ • S + Bl ■ S, the tvha ■ S terms cancel. However, trh is generated since in 
this case the site-parity, Vs, is broken. Hence we arrive at the important conclusion that a 
periodic chain with a decoupled impurity is not a stable fixed point, since relevant operators 
always will be present. 

The situation is different for the open chain. In this case the boundary operator formalism 
identifies left and right-moving operators and the chain-end spins become: 

S± oc J±. (2.3) 

Here + and — refer to the two sides of the break in the chain and J = J L = J R . Since h 
doesn't appear as a boundary operator, all perturbations of the open chain are, at most, 
marginal. A weak coupling of the two sides of the break is irrelevant. A weak coupling to an 
external spin generates (J+ + J_) -S, which is analogous to a Kondo coupling. It is marginally 
relevant for antiferromagnetic coupling and marginally irrelevant for ferromagnetic coupling. 
Hence the open chain with no decoupled impurity or with an impurity whose coupling flows 
to zero from the ferromagnetic side are stable fixed points. 

Let us now consider the stable fixed points for open spin-chains with link-parity symmet- 
ric couplings to an external S — 1 impurity, a class of models which includes the integrable 
one. The case of a simple Heisenberg coupling: 

tfhnp = J(Bx + ° L )-S (2.4) 

was discussed ear fieri. If J < 0, it renormalizes to zero leaving the open chain with a 
decoupled S — 1 impurity. If J > it is marginally relevant and we assume that it 
renormalizes to oo. This produces an open chain with two sites removed and no leftover 
impurity in the low-energy theory. 

The excitation spectrum of a long chain of length L contains towers of states with spacings 
of 0(1/ L) up to higher order corrections. This low-energy spectrum, which is a universal 
property of the fixed point, is reviewed in Ref. || for periodic and open chains. Some of the 
first few states are given in Table | of the present paper. Note that spin chains of even length, 
L, with periodic or open boundary conditions, have parity even (odd) ground-states for L/2 
even (odd). To 0(1/ L), the spectra are identical for L/2 even or odd apart from a parity 
flip for all states. Thus we see that in specifying the various fixed points we must be careful 
to specify the ground-state parity. The infinite antiferromagnetic J fixed point referred to 
in the previous paragraph corresponds to a ground-state with reversed parity compared to 
J = 0, since two spins have been effectively removed from the chain to screen the impurity. 
We will take the original chain length to have L/2 even. Thus this screened fixed point has 
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a ground-state with = 0~, where St is the spin of the state and P its parity. We label 
the renormalization group fixed point corresponding to an open chain with no parity flip as 
open + ;we label the fixed point with the parity flip as open - . The open chain with no parity 
flip and a leftover decoupled S = 1 impurity is labeled open + x (S = 1). Thus a negative J 
renormalizes to the open + x (5 — 1) fixed point and a positive J renormalizes to the open - 
fixed point. 

Let us now consider the integrable impurity model, of Eq. ( |1 . 1| ) . The Bethe ansatz 
results of Ref. ^ indicate that the magnetization has Curie form as T — > with magnitude 
corresponding to a decoupled S = 1/2 impurity. It is quite easy to see that such a fixed point 
will not arise in a link-parity invariant model from the type of analysis used above, where 
it is assumed that all couplings renormalize to oo or 0. If a single chain-spin coupled most 
strongly to the impurity it could partially screen it, leaving an effective S = 1/2 impurity. 
However, with link-parity, this type of analysis always produces changes in the effective 
spin by integer units. The other possibility is that the S = 1 impurity effectively "splits 
in half", donating an extra S = 1/2 to the chain and leaving behind a decoupled effective 
5=1/2 impurity. [See Fig. [I].] The coupling of the rest of the chain to the "donated" 
S = 1/2 must "heal", ie. renormalize to the same value as in the rest of the chain. Such a 
healing phenomena was shown to occur for an S — 1/2 impurity coupled symmetrically to 
two neighbouring sites in a chain!. However, this does not correspond to a stable fixed point 
in the present case because of the decoupled S — 1/2 impurity. A residual coupling of this 
to the healed chain is relevant, as discussed above. It generates the operator, —X'trha ■ S e s, 
where S e s is the effective S = 1/2 impurity. 

We now propose a resolution of this dilemma. Due to the very particular nature of the 
integrable Hamiltonian, Eq. ( |1 . 1|) , the relevant coupling, A', referred to above, vanishes. Of 
course, if we were to make an infinitesimal change in any of the lattice coupling constants 
near the impurity, we should expect that this relevant coupling in the fixed point Hamiltonian 
would generally become non-zero. In the next section we explore, using both RG and finite- 
size Bethe ansatz analysis, this hypothesis about the integrable impurity model itself. In 
Section IV we use RG and the modified Lanczos method to study perturbations of the 
integrable impurity model. 



III. THE INTEGRABLE IMPURITY MODEL 

We now want to study the integrable impurity model in more detail. As stated above, 
we hypothesize that it renormalizes to the unstable fixed point corresponding to the S = 1 
impurity breaking up into two S = 1/2 spins, one of which is adsorbed into the chain and 
the other of which decouples. [See Fig. [I].] We assume that the chain originally had length 
L; hence, after adsorbing the extra 5 = 1/2 variable, it has the effective length I = L+ 1. 
For this decoupling to occur, the relevant coupling of the extra spin to the chain, discussed 
in the previous section, must be "fine-tuned" to zero. The next most important coupling to 
consider is then the marginal coupling of the impurity to the periodic chain. This can be 
written: 

5H = -X(l)v(J L + J R )(0) ■ S eS , (3.1) 
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where Jl,r{x) is the spin-density of left (right)-movers and S e s is the effective impurity spin 
assumed to have size 1/2. v, the spin-wave velocity, plays the role of the velocity of light 
in the conformal field theory. Its value, v = tt/2, is known exactly from the Bethe ansatz. 
A positive A corresponds to a ferromagnetic coupling. The total spin of the left-movers is 
given by: 



S L = (1/2tt) J dxJ L} (3.2) 
and similarly for the right-movers. The total spin of the periodic chain (not including S e g) 



is 



Sdiain — Sl + Sr. (3.3) 

The /3-function for A is calculated in Appendix B of Ref. [l]. [See Eq. (BIO). There a 
positive A is antiferromagnetic] The calculation is identical to that for the Kondo problem 
and the result is: 

d\/dhxl = -A 2 . (3.4) 

Although, in that appendix we only have the left-moving part of J, the two parts of the 
interaction renormalize separately so we get the same /3-function. Note that a ferromagnetic 
coupling is marginally irrelevant. Solving, we obtain the effective coupling constant at scale 
I in terms of the effective coupling constant at scale Iq as: 

A(/) * 1 + A(^W (3 ' 5) 

Thus, if the integrable impurity model is to renormalize to the proposed fixed point, the 
marginal coupling, A must be ferromagnetic (ie. A > 0) in addition to the relevant coupling 
vanishing. There is no particular reason for the marginally irrelevant coupling, A to be 
strictly zero, and indeed we shall see that it is not. Such a marginally irrelevant coupling 
leads to corrections to the asymptotic behaviour which only vanish as 1/ln/. Consequently, 
it becomes difficult to conclude very much about the critical behaviour from finite-size 
calculations unless exponentially long chains can be studied. Fortunately, this is possible 
using the Bethe ansatz. 

A similar difficulty was already encountered for the periodic S = 1/2 chain, without any 
impurity. In that system there is a "bulk" marginal operator: 

5H = -g(l)v(8-K 2 /VS) J dxf L (x) ■ J R (x). (3.6) 

(Recall that dimension two bulk operators are marginal, but dimension one boundary op- 
erators are marginal; the difference arises from the integral over dx in the former case.) In 
this case, the renormalized coupling is given by0: 

g{l) = g(l )/(l + Ang(l ) \n(l/l )/V3). (3.7) 
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As first pointed out by CardyEl, and applied to the study of periodic Heisenberg chains in 
Ref. [13|, the effect of such marginally irrelevant couplings on the finite-size spectrum can be 
calculated in perturbation theory in the effective coupling constant. 

At a conformally invariant fixed point, excitation energies, take the formEl: 

E n -E = ^x n , (3.8) 

where Eq is the ground-state energy and the scaling dimensions, universal. Thex„'s 

for the lowest energy state of given Sl, Sr can be written: 

x n = (S L ) 2 + (S R ) 2 . (3.9) 



[See for example Ref. [13|.] The total spin multiplets are determined by the usual angular 
momentum addition rules, \Sl ~ Sr\ < 5 cna in < Sl + Sr. For an even length chain Sl and 
Sr are either both integer or both half-integer; for an odd length chain one of them is integer 
and one is half-integer. 

The ground-state energy takes the form: 

E = e l-^, (3.10) 

where eo = In 2, the ground-state energy density, is non-universal and the universal l/l 
correction is proportional to c, the conformal anomaly parameter; c = 1 for the 5 = 1/2 
chain. 

The excitation energies receive corrections in first order perturbation theory in the bulk 
marginal coupling constant, gElll!: 

Sx n = —-j=g(l)S L - Sr. (3-11) 



Note that this dot product can be determined from Sl, Sr and 5, 



chain • 



S L -S R = (1/2) [(Sl + S R ) 2 -Sl- S 2 R ] 

= (l/2)[5 chain (5 chain + 1) - 5 L (5 L + 1) - Sr(S r + 1)]. (3.12) 

The ground-state energy only obtains a correction of third order in g(lpM: 

*c=^. (3,3) 

The integrable impurity model has two marginal coupling constants, g(l) and X(l), pro- 
ducing two sources of logarithmically slow finite-size behaviour. The corrections due to g 
will be the same as for the periodic chain, given above. We now calculate the corrections 
due to the marginal boundary coupling constant, A. 

There is a first order correction to the excitation energies. Since, for A = 0, the chain is 
translationally invariant, for any eigenstate, 

< n\(J L + J R )(0)\n >= [2tt//] < n\S chain \n > . (3.14) 
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The marginal coupling of the impurity to the periodic chain, Eq. ( |3.1|) , will then give rise to 
a finite size correction of the scaling dimension of the following form 

5x n = -A(/)S c hain • S cS . (3.15) 

This can be expressed in terms of the observable total spin of the state, St where St = 
^chain + See, the spin of the periodic chain, S'chain, and the impurity spin, S e fr = 1/2, giving: 

S'chain • S ofr = (1/2)[S T (S T + 1) - S c hain(S c hain + 1) - 3/4]. (3.16) 

Combining the various terms we obtain 



E n - E 1/2 = (2nv/l) 



> 47T - 
$L + Sr ~ ^=9\1)Sl " Sr ~ -MOSchain ' ^eff 



(3.17) 



where the effective length, I, for the integrable impurity model is I = L + 1. 

The correction to the energy of the ground-state, which has Sl = Sr = S c hain = 0, 
St = 1/2, vanishes to first order in A, so let us consider next order. We can express the 
second order correction in terms of the expansion of the partition function in powers of A in 
the zero-temperature limit. This gives: 

/oo 
dr x dr 2 < 5H(t x )5H(t 2 ) > +...] (3.18) 
-oo 

Thus the correction to the ground-state energy is: 

/°° -» -> ->-,-> -> 

dr < 0\(J L + J R )(0,0) ■ S cff (0)(J L + J«)(r,0) ■ S cff (r)0 > . (3.19) 
-oo 

In lowest order perturbation theory: 

< 0|S c V0)St ff (T)|0 >= (l/3)^S c 2 ff = {1/3)3*8(8 + 1) = (1/4)^. (3.20) 
We also need: 

<0\J a L (r)J b L (0)\0>=5 ab /2(vr) 2 , (3.21) 

and the same for Jr. (Left and right are uncorrelated. See Ref. |l] for the normalization.) 
This is just the free fermion current Green's function. This is the result for an infinite 
system. To get the result for a finite system we make a conformal transformation to map 
the infinite plane onto the cylinder of circumference I, or else just work out explicitly the 
free fermion current Green's function with appropriate boundary conditions. The result is: 

< 0| J£(t)J£(0)|0 >= S ab . 1 — . (3.22) 

2[(//7rJ smh(7rt>r//)J z 

Note that the r integral in Eq. ([3 .18 ) diverges as r — ► 0. This is an ultraviolet divergence 
which would be cut off by the lattice spacing of the spin chain. It is simplest just to put 
in a cut off on the r integral, |r| > r . To evaluate the integral we change variables to 
u = tanh(7rt>r/Z), giving: 
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,)E {I = -(3/ !:)A 2 t»y / du/u 2 , (3.23) 



"0 



where w ~ ttvt /1. (We assume Z >> r .) Thus 

^ = -(3/4)A 2 i;y(l/ M o-l) 

= -3A 2 /4r + 3A 2 w(tt/4/). (3.24) 

Note that the first, cut-off dependent term, is a non-universal contribution to the l- 
independent part of Eq. The second is a correction to c: 

<5c=-^A 2 (Z). (3.25) 

The universal ground-state energy correction is second order in A but third order in g; we 
find that the A correction is much larger. Assembling the various terms the ground-state 
energy takes the form: 

^ 2 -oL- £l = -^[l + M-^(0l, (3.26) 



x/3 2' 



where we should use I — L + 1 for the integrable impurity model. Here Eq = In 2, and from 
the work of Schlottmanni e 1 = (l/2)(^(3/4) - ^(5/4)) ~ -0.4292036733 where ^ is the 
digamma function. 

The ground-state, with St = 1/2, occurs for L odd corresponding to the effective length, 
I = L + 1, being even. All states occurring for even L are regarded as excited states. 
In applying Eq. (|3.17|) , the values of g{l) and A(/) are expected to interpolate smoothly 



between even and odd I. In table (|TV|), we give all the relevant quantum numbers for the 
first few lowest energy states. Note that the states with half-integer S^am all come in pairs 
of opposite parity, obtained by interchanging the quantum numbers Sl and Sr. These pairs 
are degenerate including 0(A) corrections. Presumably they are, for more general models, 
split by corrections of higher order in irrelevant operators. However, as we shall see, for the 
integrable impurity model they remain exactly degenerate. 

In what follows, we test these formulas in two different ways. One way is to confirm that 
all energy levels that we consider are given by these formulas with the same values for g{l) 
and A(/) for a given length I = L + 1. We expect small discrepancies to occur because of 
corrections in higher orders of perturbation theory; however, these should become smaller at 
larger I. Secondly, the functions g(l) and A(Z) should be given by the lowest-order /3-function 
results, Eqs. ( |3.5| ) and ( |3.7| ) for sufficiently large /. 

Following Andrei and Johannesson@ the Bethe ansatz equations, for the integrable Hamil- 
tonian of Eq. (|1.1[), are 

' k k + il/2\ L (k k + i\ *l A fc - A; + i 

11 A A. „•> K 6 -*') 



\A k - il/2 J \A k -iJ i=i A k — A t — i 

where L is the number of S = 1/2's. The number of roots, M determines the total S z 
component through the relation S z = L/2 + 1 — M. In terms of the solutions, A^, to the 
Bethe ansatz equations, Eq. ( f3.27| ), the energy is given by 
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M 1/2 



k- 



1, Af, + (1/2)= 



One should note that the energy, Eh, of the Hamiltonian Eq. (1.1) is related to E by 
E H = E + L/A + 2/9. 

The Bethe ansatz equations are solved numerically by first making the assumption that 
the solutions occur in strings of length n. 

A"> a = A" + i(n + l - 2a)/2, a = 1, . . . , n, (3.29) 

and then solving Eq. ( |3.27| ) for the centers of the strings, A™. If v n is the number of strings 
of length n then we must have 

M = Y,nu n . (3.30) 

n 

The string assumption is then relaxed and the full Bethe ansatz equations are solved by a 
Newton- Raphson method using the string solution as the starting point. 

For chains of even length, L+l, corresponding to an odd number, L, of sites with S — 1/2 
and one spin S — 1, we determine the three lowest lying levels, Ey 2 ,E 3 / 2 and Ey 2 . Here 
the index refers to the total spin, St, of the state. The ground-state E 1 / 2 corresponds to 
a solution with (L + l)/2 real roots and no strings, E 3 / 2 is a solution with (L — l)/2 real 
roots and Ey 2 has (L — 3)/2 real roots plus a 2-string at x ± i(jr/2 + 5) where 5 is a small 
positive number quickly approaching zero for long chains, and x is non-zero. The results for 
these three levels are given in Table ||. For chains of odd length, L + l, corresponding to 
an even number, L, of sites with 5=1/2 plus one site with spin S = 1, we determine the 
two levels E\ and E 2 . E\, the lowest lying state, is a solution of the Bethe ansatz equations 
with L/2 real roots. We were not able to obtain the first excited state, Eq, by the Bethe 
ansatz scheme. The next excited level, E 2 , has L/2 — 1 real roots. The results for E\ and 



E 2 are summarized in Table Q l 



These states exhibit some remarkable degeneracies, for finite L. Eq. ( |3.27| ) is invariant 
under the operation, {A^} — > {— A^}. Thus, in cases where the set of roots is not symmetric 
about 0, a pair of degenerate solutions is obtained, if we assume that the corresponding 
wave-functions are linearly independent as can be verified for short chains. This is the case 
for all solutions discussed above with even L and also for (1/2)*. The even L degeneracies 
can be understood from the picture of the RG fixed point corresponding to a periodic chain 
of (L + 1) S — 1/2's and a decoupled S = 1/2. The periodic odd-length chain has a large 
exact degeneracy for finite L. This follows from the fact that the ground-state does not have 
zero momentum. Instead it consists of two degenerate doublets, (l/2) ± with momentum 
±/co- As L —>■ oo, we expect ko — ► n/2. By forming linear combinations of these states 
we can form positive and negative parity eigenstates. All low-lying states have momentum 
close to ±7r/2 and consequently also come in parity doublets. In the conformal field theory 
picture, the parity doubling arises from the fact that the periodic chain of odd length, L + l 
has Sl integer and Sr half-integer or vice versa. The corrections to the excitation energies 
of first order in g and A, Eq. ([3.171) do not lift the degeneracy. This must be true to all 



orders in g and all other irrelevant bulk operators for the periodic chain. On the other hand, 
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we expect that higher order corrections in irrelevant boundary operators will, in general, lift 
the degeneracy since momentum is not, in general, a well-defined quantum number for the 
impurity system. Remarkably, this does not happen for the integrable impurity model and 
the degeneracy remains exact. Even more surprising is the degeneracy for odd L. In this case 
the corresponding periodic chain has even length, L + 1 and does not exhibit any exact finite 
length degeneracies. Nonetheless, such a degeneracy occurs for the (l/2) ± * states. We sec 
from Table [TV] that this occurs via a cancellation of the 0(g) and 0(A) marginal corrections. 
At large I, using \(l) — > l/ln(/), g(l) — > y/3/Air ln(7) the correction to x n becomes: 

Sx n = -[S L ■ S R + ,S chain • S cS ]/M. (3.31) 

[S L -S R + S chauin -S eS ] has the value -3/4 for both the (1/2)*+ and (l/2)*~ state. The splitting 
of these levels by the marginal bulk operator is cancelled by the splitting due to the marginal 
boundary operator! Furthermore, this cancellation must be exact to all orders in irrelevant 
bulk and boundary operators. This suggests the existence of some sort of hidden symmetry 
in the integrable impurity model reminiscent of the Yangian symmetry discussed recently 
for the 1/r 2 periodic Heisenberg chainlll. 

We can now try to extract the coupling constant, X(l) defined in Eq. ( |3.1| ), as determined 
from the five levels, Eyz, E 3 / 2 , E*T 2 , E\ and E 2 . The bulk marginal coupling constant, g(l), 
has already been determined from finite-size analysis of long periodic chains without the 
impuritylll. We shall use this as input and determine the boundary marginal coupling, A(/) 
from our data on the integrable impurity model. Thus the only free parameter is A(Z). 
Note that (1/2) is the ground-state so we fit its energy to Eq. ( |3.26| ). The other four are 
excited states so we fit their excitation energies by Eq. (|3.17f ). The different estimates of 
X(l) obtained from the different energy levels are shown in Fig. [| Notice, first of all that 
A is positive, corresponding to a ferromagnetic, marginally irrelevant coupling of S e g to the 
chain. As expected the different estimates of A(Z) are all approximately the same. As I 
increases, the estimated couplings approach each other and get smaller. 

We compare A(/) to the one loop (3 function result Eq. ( |3.5| ), A rg , in Fig. [3|. It is seen 
that the first-order /3-function result is valid for chains longer than I ~ 100, indicating that 
our perturbative results should be meaningful for chains of this length or longer. A similar 
plot for the other coupling g is shown in Fig. |j. Again we see that the one loop (3 function 
gives good results for chains longer than a few hundred. As shown in Table [TV] the excited 
state, (l/2)* + , does not receive a correction to its excitation energy of first order in A. We 



can therefore use it to extract an estimate of g(l) independent of the results from Ref. [13. 
This estimate , f?i/2 (0 , is also shown in Fig. f|. It compares nicely to the f3 function results 
as well as to the results obtained for the pure periodic S — 1/2 chain in Ref. |13[ [The 



corresponding plot in Ref. [13| contained a numerical error in the calculation of g rg .} 

From the above results we conclude that the spectrum of the integrable impurity model 
indeed asymptotically becomes that of a periodic chain with L + l spin-l/2's and a decoupled 
S c s = 1/2. The impurity spin S = 1 has effectively split in half. The spectrum corresponds 
to column 1 in Table |. 
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IV. GENERAL HAMILTONIANS 



In this section, we consider the effects of perturbing the integrable impurity Hamiltonian 
Hi nt , Eq. ( p..l|) , away from its integrable form. A general phase diagram is hypothesized 
and is tested using the modified Lanczos method for general Hamiltonians on chains with 
L < 20. To be concrete, we consider the following two-parameter set of Hamiltonians: 

1 L-1 

H = - ^2 Bi ■ a i+1 + Hi 
4 i=i 

2 7 1 7 7 

Hi = +a L )-S+ -{a x -S,a L -S}- -a t ■ a L ] + Mb x + a L ) ■ 5. (4.1) 

y z o z 

Here 5 is the 5 = 1 impurity. J\ = 1, J 2 — is the integrable impurity model, and J\ = 
corresponds to the models discussed in Ref. []. J\ = J 2 = is the simple open chain with 
a decoupled 5=1 impurity. By combining our understanding of the phase diagram near 
\Ji\ = 0, \Ji\ — > oo and the integrable impurity model, we hypothesize a general phase 
diagram. 

The vicinity of |Jj| = can be readily analyzed. Near the open chain fixed point the 
terms {B\ ■ S,ctl • 5} and B\ ■ ai in Eq. (|4.1|) contain a product of both chain-end spins 
<7i and o l- In the continuum limit they correspond to a dimension 2 irrelevant operator 
containing the product of J + and J_. Thus for small Jj we may approximate the impurity 
part of the Hamiltonian, Hj, by (2Ji/9 + J2/2)(a + a^) ■ 5. This is marginally irrelevant 
for (J 2 + 4Ji/9) < 0. Thus, for (J 2 + 4Ji/9) < 0, the system renormalizes to the fixed 
point consisting of the open chain with a decoupled S = 1 impurity which we denote by 
open + x (S = 1). For (J 2 + 4Ji/9) > the coupling is antiferromagnetic and therefore 
marginally relevant. We expect it to renormalize to oo, screening the impurity and effectively 
removing two sites from the chain. This produces a parity flip, giving the open - fixed point. 
(Recall from Sec. II that the removal of two spins from the open chain by the screening 
process reverses the parity.) This transition is oo— order since it is driven by a marginal 
operator, i.e. the cross-over length scale diverges exponentially as (2Ji/9 + J%/2) — > from 
positive values. 

Let us next consider what happens in the limit | J$| — ► oo. In this limit we must find the 
ground-state of the three-spin system, S, <?i, c?l, of Hi. Depending on the Jj's, the ground- 
state has spin and parity + , 1 + or 2 + . These wave-functions are depicted schematically 
in Fig. H as are the regions of stability of the three states. The + case leads directly to 
a stable fixed point when we include the relatively weak coupling of these three spins to 
the rest of the chain. These couplings are irrelevant, leaving the open chain fixed point 
open - . The 2 + impurity state, however, is unstable because the effective 5 = 2 impurity is 
coupled antiferromagnetically to the chain. Assuming this coupling flows to oo, the impurity 
is partially screened, leaving an effective 5 = 1 impurity with an a ferromagnetic coupling 
to the rest of the chain. Hence, it will renormalize to 0. Since four sites are involved in 
producing the effective 5=1 decoupled impurity, there is no parity flip. Thus we obtain 
the open + x (5 = 1) phase. Of special interest is the 1 + ground-state for the three-spin 
complex. d\ and <7l form an 5 = 1 spin, symmetric under interchanging these two spins. 
This then couples to the 5=1 impurity to form a state of total spin 1. Note that neither 
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the minimal spin (0) nor the maximal spin (2) occurs. The biquadratic term in iff, involving 
the anticommutator, is necessary to assure that this happens. Thus we begin to see how 
the integrable impurity model can have different behaviour than the simpler one discussed 
above, with J x = 0. The 1 + state of the three-spin cluster does not correspond to a stable 
fixed point since the rest of the chain is coupled antiferromagnetically to the effective 5 = 1 
impurity, and such a coupling is marginally relevant. We expect this coupling to the next 
two spins in the chain, 02 and Bt-x t° a ls° renormalize to 00, screening the effective impurity 
and leading to an open chain fixed point with no leftover impurity spin. In this case, four 
chain-spins are involved in screening the impurity in a parity-symmetric way, so that the 
stable fixed point is open + . 

For large | Jj|, we expect the above three different stable phases to occur, with the phase 
boundaries asymptotically approaching those of the three-spin system, as drawn in Fig. |6|. 
Note that the open + x (5 = 1) phase is equivalent to Jj = 0, an open chain with a decoupled 
5 = 1 impurity; we may think of this entire phase as renormalizing to the origin. The open + 
and open - phases, however, are characterized by impurity couplings renormalizing to 00. 

We have established, in the previous section, that the integrable impurity model renor- 
malizes to a fixed point in which the 5 = 1 impurity effectively splits in half, one extra 
5 = 1/2 being absorbed by the chain and the other decoupling. This can only occur if 
no relevant operator connects the periodic chain to the decoupled 5 = 1/2 spin. We now 
analyze the effect of small perturbations around the integrable impurity model by consider- 
ing the periodic chain with both marginal and relevant couplings to the effective 5=1/2 
impurity, i.e. we consider the continuum limit Hamiltonian consisting of the k = 1 WZW 
model with two local perturbations: 

SH = -[X'tvha + X(J L + J R )\ ■ S cS . (4.2) 

The integrable impurity model has A' = 0, A > 0; ie. the relevant coupling vanishes and the 
marginal one has the irrelevant sign. An infinitesimal perturbation of the integrable impurity 
model will, in general, produce a non-zero A . The resulting behaviour was discussed in 
Ref. ^, in the context of a simple coupling of an 5 = 1/2 impurity to a single site in a periodic 
chain. We expect that A will renormalize to ±00 beginning from an infinitesimal positive (or 
negative) value. The negative case, corresponding to an antiferromagnetic coupling, leads 
to screening of 5 e g by a single site in the chain. The stable fixed point is an open chain 
with one site removed, open" 1 ". For A > 0, the ferromagnetic case, S e g and the site to which 
it is coupled form an effective 5 = 1 impurity. However, this is not a stable fixed point. 
The 5=1 effective impurity is coupled antiferromagnetically to two neighbouring spins. 
We expect this coupling to renormalize to 00, screening the effective impurity. Once again 
the stable fixed point is an open chain. However, in this case, three chain spins are involved 
in the screening process and get removed from the effective open chain at the stable fixed 
point. As discussed in Section II, the parity of all low-energy states is flipped relative to 
the case where a single chain-spin is removed. Thus the stable fixed point in this case is the 
open - . We see that the unstable critical point, A = A = 0, to which the integrable impurity 
model renormalizes separates the stable open" 1 " and open - phases. 

We now turn to a discussion of the order of the phase transition separating the open + 
and open - phases. We expect that a second order critical line will exist for a finite range 
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of positive A with A' = governed by the A = A' = critical point. On the other hand, if 
the marginal coupling A < 0, then it is relevant and renormalizes to large values. In this 
case, the simplest assumption is that A — > — oo; otherwise we would be forced to postulate 
another non-trivial critical point. In general, when impurity couplings renormalize to oo we 
expect a first-order phase transition. The reason is that we can then ignore any couplings 
of the impurity complex to the rest of the chain. In this particular case we can consider 
only the 5 = 1/2 impurity and three chain spins. This cluster of four S = 1/2 's has a + 
or - ground-state depending on the various couplings. The phase transition in this limit 
is a simple level-crossing in the four-spin system and is therefore first order. The critical 
point, Pi, at A = A = is on the phase boundary between open" 1 " and open - phases and 
separates the second from first-order transition lines. The integrable impurity model, which 
was shown in the last section to have a non-zero positive A, lies on the second-order part of 
the phase boundary. 

Now we attempt to combine our information about the small Jj, large Jj regions and 
the vicinity of the integrable point. Our large J« analysis tells us that there are three 
stable phases, open - , open + and open + x (5 = 1). Our analysis of the vicinity of the 
integrable impurity model tells us that it should be on the open - — open + phase boundary. 
It renormalizes to a critical point, Pi in Fig. || where this transition changes from first to 
second order. 

The open - phase occurs when we increase J 2 from Pi. This is to be expected because J 2 
corresponds to a coupling of S e s to the two nearest neighbours, <?i, c?l, not to the adjacent 
spin c?o, as shown in Fig. [?|. Antiferromagnetic J2, i.e. A' > 0, leads to a screening of S c s 
by c? , a 1 and a^. The removal of c?o, d\ and &l from the open chain leaves a chain with 
L — 2 sites and thereby implies a parity flip. In Fig. |3| we show a Lanczos calculation of some 
low-lying states for length 20, as a function of J 2- The open - spectrum shown in Table | 
occurs for positive The Lanczos results are discussed in more detail below. 

We see that there must be another multi-critical point in the phase diagram where all 
three stable phases meet, P2 in Fig. [| This point is presumably not at the origin since 
there is only one marginal operator in the vicinity of the origin, (2^/9 + J 2 /2)(cx + &i) • S, 
as discussed above, so we only expect two phases to meet at that point. We hypothesize 
that this higher multi-critical point corresponds to an open chain with two S = 1/2 effective 
impurities, Si and S2, decoupled from the chain and from each other, ie. the original S = 1 
impurity effectively breaks up into two S — 1/2 impurities, with everything decoupled at 
the multi-critical point. There is one relevant coupling and two marginal ones at this critical 
point. We write these, in the lattice model as: 

5H = eSi ■ S 2 + eiSi ■ {B\ + a L ) + e 2 S 2 • (c?i + a L ). (4.3) 

[See Fig. |9|.] We now analyze the phase diagram for this model. Since the impurity spins 
have zero scaling dimension, the coupling e is highly relevant. Assuming that a non-zero e 
renormalizes to ±00, the two impurities lock into a singlet leaving the open chain fixed point, 
open + , for e > 0. For e < they lock into an effective S — 1 impurity. If both couplings 
ei < (ferromagnetic), this S = 1 impurity decouples, giving the open + x (S — 1) phase. 
On the other hand if at least one of the couplings is antiferromagnetic, then the effective 
S = 1 impurity is screened, giving the open - phase. The open" 1 " x (S = 1) — open - phase 
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transition is equivalent to the one discussed in the second paragraph of this section. It is 
therefore of oo-order. Note that the open + x (S 1 = 1) to open" 1 " phase transition is controlled 
by the multi-critical point, P2 at e = = 0. This transition is governed by the relevant 
coupling constant e, which only involves the two impurities and not the rest of the chain. 
Thus it corresponds to a simple level-crossing in this two-spin system and so should be first 
order. We hypothesize that, as we move along the critical line, where e = 0, one of the 
marginal couplings, say ex, changes sign at the multi-critical point, while the other remains 
ferromagnetic. The transition between open - and open" 1 " phases, for e\ > 0, is governed by 
the behaviour at e\ of 0(1). We know that a single S = 1/2 impurity with such a coupling 
to an open chain will get absorbed by the chain, ie. the defect heals and the fixed point is the 
periodic chain with an extra spin. Thus it is plausible that this phase boundary is second- 
order and renormalizes to a periodic chain with a single decoupled S = 1/2 impurity. Note 
that the continuation of the open" 1 " x (S = 1) — open" 1 " phase boundary is the open - — open 4 " 
phase boundary, but the order of the transition changes from first to second, as drawn in 
Fig. ^. This follows from the fact that both phase boundaries are governed by the vanishing 
of the relevant coupling constant, e. The system renormalizes to Pi by one of the impurities 
being absorbed into the chain. It renormalizes to the origin by the two impurities locking 
into a decoupled S — 1 impurity. 

Altogether there are five different critical points; three occurring at finite coupling, P 1; P 2; 
open + x (S = 1) and two at infinity, open 4 ", open - , as shown in Fig. || Various sections of the 
transition lines are first, second or infinite order. The detailed shape of the phase boundaries 
depicted schematically in Fig. ^| is not known. What is known is (i) the asymptotic slope 
of the three phase boundaries at |Jj| — > 00, (ii) the slope of the phase boundary at the 
origin and (iii) the fact that the integrable impurity model at J\ — 1, Ji = lies on the 
open - — open 4 " phase boundary. 

We now discuss our numerical results on chains of length L < 20. We emphasize at 
the outset that we are fighting finite-size corrections that vanish logarithmically slowly from 
two sources: the bulk marginal coupling, g of Section III which is present everywhere in the 
phase diagram and the marginal boundary operator A which is present in some parts of the 
phase diagram. Thus we can only expect our results to be of anecdotal value. 

In Fig. we present the scaled energy gaps, L(E — E 1 +)/2ttv, for the four states with 
quantum number O 4 " and l 41 as a function of J\ with J2 fixed at J2 = 0. In Fig. [It] the 
integrable impurity model thus corresponds to J\ = 1. 

First, let us consider what happens as we increase J\ away from its value, J\ — 1 at 
the integrable impurity model. We see that the + — l - gap drops rapidly with increasing 
J\. This is to be expected since, according to Fig. [J] the system is in the open 4 " phase. 
Note however, that the + state is not the ground-state, even for J\ « 4, for L < 20. We 
do expect that it would become the ground-state for sufficiently large L, for any J\ > 1. 
In Fig. [10] we show results for two different chain lengths L = 8 and L = 20. As can be 
seen there is some evidence that the l - — + gap indeed is closing with increasing L. The 
scaled gaps between the O 4 " — l 4 " states become asymptotically degenerate at J\ = 1, the 
integrable impurity model. As seen in Fig. ^ the + level has the most negative slope as 
a function of J\. It is then plausible that as the gap closes at J% = 1, the + state cross 
the l - state for J\ > 1. Even for large J\ this level-crossing only takes place at large L. 
The reason is that as J\ — ► 00 for fixed L, we obtain the ground-state of the three-spin 
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cluster, 1 + , times the ground-state for the rest of the (open) chain with two sites removed 
resulting in the 1~ ground-state. Eventually, if J\ is kept fixed, the effective S = 1 impurity 
is screened for long enough chains, producing the + ground-state, but this process proceeds 
logarithmically slowly. Finite-size scaling analysis, although not very reliable due to the 
marginal operators, seems to indicate that J 2 is relevant, as expected. 

Now consider what happens as we decrease J\. At J\ = we obtain, approximately the 
open + x (S — 1) spectrum, shown in Table |, as expected. Note also, the minimum in the 
1 + — - gap which occurs for Ji « .7. With increasing L we expect a crossing of these two 
levels to occur. Asymptotically the - state should lie below the 1 + state for all Ji such 
that < J 1 < 1, since this region should be in the open - phase. 

We now consider the effect of varying J 2 away from with J\ held fixed at its integrable 
value, J\ = 1; [see Fig. The open" spectrum of Table | is obtained for large positive 
J 2 , as expected. Although we expect to be in the open + x (S = 1) phase for sufficiently 
large J 2 < 0, this is not obvious from Fig. The problem is that, for J 2 — > —00 for fixed 
L the three-spin complex has an S = 2 ground-state giving a 2~ ground-state including the 
decoupled open chain. The screening of this 5 = 2 effective impurity down to S — 1 is 
logarithmically slow. Indeed, the open + x (S = 1) spectrum is best approximated in this 
region for J 2 ~ — -3. This may correspond to crossing the open - — open + x (S = 1) phase 
boundary, in Fig. ^ where the marginal boundary coupling vanishes. 

According to the "(/-theorem" the "ground-state degeneracy", g, decreases under renor- 
malization from less stable to more stable fixed points.H The value of g for a periodic chain§ 
is 1 and for an open chain is 1/ y/2. These values must be multiplied by the degeneracy of 
the decoupled impurity at each critical point. Thus, at P 2 , g = 4/v2; at Pi, g = 2, at the 
origin, g = 3/^/2 and at the open" 1 " and open - critical points at 00, g = 1/V2. We see that, 
in all cases, the (/-theorem is obeyed. A pictorial summary of the five fixed points and the 
corresponding values of g are given in Fig. [II]. 

V. CONCLUSIONS 

The apparent contradiction between the integrable impurity model and our RG analysis 
is explained by the fact that this model renormalizes to an unstable critical point correspond- 
ing to a periodic chain with a decoupled S — 1/2 impurity. This has been shown rather 
convincingly for the integrable impurity model itself from finite-size analysis of chains of 
length L < 5, 000. A general phase diagram has been proposed and supported by finite-size 
analysis of general models on chains with L < 20. 
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TABLES 



TABLE I. The low temperature spectra showing only the four states + ,0 ,1 + ,1 for the 
various fixed points. L is divisible by 4, x = L(E n — Eq)/2ttv. 



X 


Integrable Model: 


Open Chain 


Open Chain 


Open Chain 




Periodic Chain 


of L — 2 sites 


of L or L — 4 sites 


of L sites with 




of L + 1 sites with 


with singlet 


with singlet 


decoupled spin-1 




decoupled spin- 1/2 impurity 






impurity 




Pi 


open - 


open" 1 " 


open + x (S = 1) 


7/4 




0+ 


o- 




5/4 




1- 


1+ 


0+ 


3/4 




1+ 


1- 


o-,i- 


1/4 


o+,o-,i+,i- 


o- 


0+ 


1+ 



TABLE II. Spectrum of the integrable impurity model for even chains, corresponding to an 
odd number, L, of sites with S=l/2 plus one spin S=l. The levels shown are the ground- state, 
#1/2, which has S=l/2, the first excited state, £"3/2, which has S=3/2, and the second excited 
state, E*, 2 , which has S=l/2. 



L 


-#1/2 


-E3/2 


771* 

-^1/2 


3 


2.666666666667 


2.000000000000 


0.833333333333 


5 


4.000000000000 


3.535183758488 


2.788431731592 


7 


5.360920843433 


5.000000000000 


4.473627449978 


9 


6.732050807569 


6.435509583680 


6.035934151803 


11 


8.108179626193 


7.855772506636 


7.536851526445 


13 


9.487154267776 


9.267023613165 


9.003326289075 


15 


10.867912224042 


10.672479361373 


10.448698007398 


17 


12.249865197376 


12.073976382835 


11.880255773559 


19 


13.632659805114 


13.472643412759 


13.302298857367 


21 


15.016070414566 


14.869214052794 


14.717519700013 


99 


69.058063860390 


69.021685474014 


68.996582070614 


199 


138.369221288792 


138.350334224364 


138.339077391494 


399 


276.996826009204 


276.987093406252 


276.981995627127 


599 


415.625640252833 


415.619054638227 


415.615834877143 


799 


554.254762484295 


554.249776039106 


554.247448415855 


999 


692.884009113116 


692.879992507487 


692.8781812698 


2047 


1419.301863049853 


1419.299863460628 




4999 


3465.472116743763 


3465.471281523977 
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T A T3T I? TTT 

lArSLij 111. 


Spectrum of the integrable impurity 


model for odd chains 


, corresponding to an 


even number, L. 


of sites with S=l/2 plus one spin S= 


=1. The levels shown 


are the ground- state, 


E\, which has S : 


=1, and the second excited state, E2, 


which has S=2. 




L 


-Eh. 




—E2 


2 


1.333333333333 




0.000000000000 


4 


2.951367322083 




1.767591879244 


6 


4.413722666901 




3.489392593998 


8 


5.839581489564 




5.090273141761 


10 


7.250399759403 




6.621632325853 


12 


8.653468913494 




8.111911768408 


14 


10.051996142025 




9.576268507776 


16 


11.447623604652 




11.023300473744 


18 


12.841282703621 




12.458202998467 


20 


14.233542424308 




13.884288605831 


100 


69.731013913337 




69.651755717692 


200 


139.052006805828 




139.011131675689 


400 


277.684686429578 




277.663763091977 


600 


416.315227987524 




416.301121655718 


800 


554.945223280210 




554.934568453894 


1000 


693.574997489704 




693.566430520955 


2046 


1418.607647701090 




1418.603401388401 



TABLE IV. 


The quantum 


numbers 


Sl, Sr, 


*S c hain> S*L " Sr, S'chain 


• S e s, and x 


for the five 


levels, describing the integrable impurity 


model. 










S L 


Sr 


'S'chain 


Sl ■ Sr 


"S'chain " 'S'cflF 


X 


L+l Odd 














2+,2- 


1 


1/2 


3/2 


1/2 


3/4 


5/4 


0+,0~ 


1/2 





1/2 





-3/4 


1/4 


1+1- 


1/2 





1/2 





1/4 


1/4 


L+l Even 














1/2*+ 


1/2 


1/2 





-3/4 





1/2 


1/2* - 


1/2 


1/2 


1 


1/4 


-1 


1/2 


3/2- 


1/2 


1/2 


1 


1/4 


1/2 


1/2 


1/2+ 
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FIGURES 



FIG. 1. The S = 1 impurity effectively "splits in half", donating an extra S = 1/2 impurity to 
the chain. 

FIG. 2. The coupling constant, Aj, for the five levels, £1/2, -£"3/2; ^1/2 > -^1 an< ^ ^2> as a function 
of the effective chain length, I = L + 1. 

FIG. 3. The average of the five coupling constants, A av , A = (\\/2 + A 3 / 2 + A^ 2 + Ai + A2)/5, 
as a function of the effective length, / = L + 1, compared to the one- loop renormalization group 
prediction, A rg = A av (/o)/(l + A av (/o) ln(///o))- In this plot we have used Iq = 1000. 



FIG. 4. The coupling gjZ and the average of the coupling, g av , from Ref. 13, as a function of 
the effective chain length, I = L + 1, compared to the one-loop renormalization group prediction, 
9rg = 9av(^o)/(l + 47T(7 av (7o) ln(l/lo) /^/3). In this plot we have used Iq = 1000. 

FIG. 5. Schematic drawing of the ground-state wave- functions for the three-spin system de- 
pending on Ji . "0" denotes the S z = state of the 5 = 1 impurity. 



FIG. 6. The conjectured phase diagram for the general Hamiltonian, Eq. fl4.l|) , in the parameter 
space Ji, J2. The couplings A, A' are defined in the vicinity of the multi-critical point Pi in Eq. ( |4.2[ ). 
The couplings e, are defined in the vicinity of the multi-critical fixed point, P2, in Eq. (|4.3|) . 



FIG. 7. Sefj couples with strength J2 to the two nearest neighbour spins, ai and ai- Antifer- 
romagnetic J2 leads to screening of S e s by three chain spins and hence to the open - phase. 

FIG. 8. Scaled energy gaps, L(E — E 1 +)/2irv, as a function of J2 (Ji = 1, L = 20). 

FIG. 9. The unstable fixed point P2 occurring in the phase diagram, Fig. ^. 

FIG. 10. Scaled energy gaps, L(E — £' 1 +)/27rf , as a function of J\ (J2 = 0, L = 8, 20). 

FIG. 11. The five fixed points and the RG flows occurring in Fig. ^. "0" denotes the S z = 
state of the S = 1 impurity. 
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